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We consider an electronic model for realizing the Spin Bose-metal (SBM) phase on a 2-leg trian- 
gular strip - a spin liquid phase found by D. N. Sheng et a/.[Phys. Rev. B 79, 205112 (2009)] in a 
spin-1/2 model with ring exchanges. The SBM can be viewed as a "C1S2" Mott insulator of elec- 
trons where the overall charge transporting mode is gapped out. We start from a two-band "C2S2" 
metal and consider extended repulsion motivated by recent ab initio derivation of electronic model 
for K-ET spin liquid material [K. Nakamura et ai, J. Phys. Soc. Jpn. 78, 083710(2009)]. Using 
weak coupling renormalization group analysis, we find that the extended interactions allow much 
wider C2S2 metallic phase than in the Hubbard model with on-site repulsion only. An eight-fermion 
Umklapp term plays a crucial role in producing a Mott insulator but can not be treated in weak 
coupling. We use Bosonization to extend the analysis to intermediate coupling and study phases 
obtained out of the C2S2 metal upon increasing overall repulsion strength, finding that the SBM 
phase is a natural outcome for extended interactions. 

PACS numbers: 



I. INTRODUCTION 

There has been much recent interest in gapless 
spin liquids stimulated by the appearance of sev- 
eral experimental candidates, including two-dimensional 
(2D) triangular lattice based organic compounds^^^ 
k-(ET) 2 Cu 2 (CN) 3 and EtMe 3 Sb[Pd(dmit) 2 ]2 and 3D 
hyper- kagome material 6 - Na^^Os- One line of theo- 
retical ideas considers states with a Fermi surface of 
fermionic spinonsi 7 -^^ For the 2D spin liquids, such a 
state arises as a good variational wavefunction 7 - for an ap- 
propriate spin model with ring exchanges; it is also an ap- 
pealing candidate for the Hubbard model near the Mott 



transition 
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Theoretical description of such states 
leads to a U(l) gauge theory (see Ref. [ID for a review). 

Variational studies are not sufficient to prove that a 
given state is realized and the gauge theory is not fully 
reliable in 2D. Driven by the need for a controlled the- 
oretical access to such phases, Ref. [13 considered the 
Heisenberg plus ring exchanges model on a two-leg tri- 
angular strip (so-called zigzag chain). Using numerical 
Density Matrix Renormalization Group (DMRG), Varia- 
tional Monte Carlo (VMC), and analytical Bosonization 
treatment, Ref. [3 found a ladder descendant of the 2D 
spin liquid in a broad range of parameters and dubbed 
this phase "Spin Bose-Metal" (SBM). The name refers to 
metal-like itinerancy present in the spin degrees of free- 
dom, while there is no electric transport to speak of in the 
spin-only model, which is bosonic model microscopically. 

A low energy field theory^ for the zigzag SBM phase 
can be obtained by employing bosonization to analyze the 
spinon-gauge theory (the slave particle approach also un- 
derlies the VMC trial states). An alternative derivation 
of the SBM theory is to consider an interacting model of 
electrons hopping on the zigzag chain and to drive a tran- 
sition from a two-band metal to a particular Mott insu- 
lator. Specifically, let us start in the metallic phase with 



two gapless charge modes and two gapless spin modes 
- so-called "C2S2" metal. We can imagine gapping out 
just the overall charge mode to obtain a "C1S2" Mott 
insulator with one gapless "charge" mode and two gap- 
less spin modes, where the former represents local current 
loop fluctuations and does not transport charge along the 
chain. This is precisely the SBM phase. If one thinks of 
a spin-only description of this Mott insulator, the gapless 
"charge" mode can be interpreted as spin singlet chiral- 
ity mode. Ref. also identified a valid Umklapp term 
that can drive the electron system to the C1S2 phase. 

In this paper, we focus on realizing such scenario for 
the SBM in explicit and realistic electronic models. Hub- 
bard model on the zigzag chain (ti — t-2 — U chain) 
has received much attention. 15 i 16 i 17 i 18 i 19 i 20 For free elec- 
trons, the two-band metal appears for £2/^1 > 0.5. How- 
ever, in the case of Hubbard interaction, weak coupling 
approac h 16 ! 17 finds that this phase is stable only over a 
narrow range £2/^1 € [0.5,0.57], while a spin gap opens 
up for larger tijt\. The Umklapp that can drive a tran- 
sition to a Mott insulator requires eight fermions and is 
strongly irrelevant at weak coupling. Prior wor k 15 i 16 i 21 
focused on the spin-gapped metal and eventual spin- 
gapped insulator for strong interaction, while the C1S2 
spin liquid phase was not anticipated. 

There have also been numerical DMRG studies of the 
Hubbard model i 17 i 18 ' 19 i 20 The focus has been on the 
prominent spin-gapped phases and, in particular, on the 
insulator that is continuously connected to the dimerized 
phase in the J± — J 2 Heisenberg model, which is appropri- 
ate in the strong interaction limit U ti,t2- The C2S2 
metallic phase and possibility of nearby spin liquid on the 
Mott insulator side in the Hubbard model have not been 
explored. We hope our work will motivate more studies 
of this interesting possibility in the Hubbard model with 
intermediate U close to the C2S2 metal. 

Since the C2S2 metallic phase is quite narrow in the 
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Hubbard model, we would like to first widen the C2S2 
region. To this end, we explore an electronic model with 
extended repulsive interactions.— Such interactions tend 
to suppress instabilities in the electronic system, sim- 
ilar to how long-ranged Coulomb repulsion suppresses 
pairing in metals. They are also more realistic than the 
on-site Hubbard, particularly for materials undergoing a 
metal-insulator transition where there is no conduction 
band screening on the insulator side. Thus, recent ab 
initio model construction for the k-(ET)2Cu2(CN)3 ma- 
terial found significant extended interactions in the cor- 
responding electronic model on the half-filled triangular 
latticed 

Applying weak coupling renormalization group (RG) 
approach to the zigzag ladder system , 15 ' 16 ' 21 ! 25 we in- 
deed find that extended interactions open a much wider 
window of the C2S2 metal phase. Building on this, we 
then use bosonization approach to explore a transition 
to a Mott insulator upon increasing the overall repulsion 
strength. We find that such longer-ranged interactions 
can drive the system into the Cf S2 spin liquid Mott in- 
sulator rather than a spin-gapped insulator. This bodes 
well for finding spin liquid phases in more realistic elec- 
tronic models for materials near the metal-insulator tran- 
sition. 

The paper is organized as follows. In Sec. [Hi we set 
up the weak coupling R G 16 ' 21 ' 25 and open a much wider 
window of the metallic C2S2 phase by introducing real- 
istically motivated longer-ranged repulsion. In Sec. IIII) 
we use bosonization to extend the analysis to intermedi- 
ate coupling. We gradually increase the overall repulsion 
strength and determine thresholds for a Mott transition 
driven by the eight-fermion Umklapp term and also for 
spin gap instabilities, thus mapping out phases neighbor- 
ing the C2S2 metal. In Sec. HVl we summarize our results 
and conclude with some discussion. 



II. WEAK COUPLING ANALYSIS OF ti t 2 
MODEL WITH EXTENDED REPULSION: 
STABILIZING C2S2 METAL 
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FIG. 1: Electron band for t 2 > 0.5ti has two occupied Fermi 
sea segments. This is free fermion C2S2 metal. 



In weak coupling, the kinetic energy Eq. (TT]) gives free 
particle dispersion 



e(k) = -2<i cos(fc) - 2t 2 cos(2fc) . (3) 



For £2/^1 > 0.5, there are two sets of Fermi points at 
wavevectors ±kpi and ±kp2 as shown in Fig. [T] We 
adopt the same conventions as in Ref. [Til Fermions near 
kpi and kp2 are moving to the right, and the correspond- 
ing group velocities are v\, V2 > 0. Electrons are at half- 
filling, which implies kp 1 + kp2 = —ir/2 mod 2tt for the 
choices as in Fig. [T] 

The electron operators are expanded in terms of con- 
tinuum fields, 



A. Setup for two-band electron system 

We consider half-filled electronic t\ — £2 model with 
extended interaction described by the Hamiltonian H = 
Hq + Hy , with 

H = - Yj [hcUx)c a (x + 1 ) + t 2 ci(x)c a {x + 2) 

+H.c] , (1) 

H v = - V(x - x')n(x)n(x') . (2) 

Here c(c^) is fermion annihilation (creation) operator, x 
is a site label on the one-dimensional (ID) chain, and 
a =|, J. is a spin index; n{x) = ^2 a c\(x)c a (x) is electron 
number on the site. 



c a {x) = Y, elPkFaXc P™ > ( 4 ) 

P,a 



with P = R/L = +/— denoting the right and left movers 
and a = 1, 2 denoting the two Fermi seas. 

Four-fermion interactions can be conveniently ex- 
pressed in terms of chiral currents; 14 ' 21 ' 25 



Jpab = Y C Paa Cpba ' ( 5 ) 
a 

Jpab = Y C Pa a ^Y- C PbP ■ ( 6 ) 

a,0 
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Most general four-fcrmion interactions can be written as, 

^■RL = X/ { W ab^RabJLab + A^ fc J R aa JLbb) , (7) 
a . b 

T~LrL = - (WabJRab ' J Lab + Ktb^ Raa 1 ^Ltb) ,(8) 

o . /.) 

^chiral = 2 ^ ^ ^aa^RaaJRaa + J LaaJ Lao) 
a 

+ Ci 2 (Jr11 JR22 + Jl\\Jl22) , (9) 
X ^ , ^ aa [J Raa ' J Raa ~t~ «7 Laa ' J Laa) 



chiral 



Ci 2 {Jrii ■ Jrti + >/lii ' */l2 



(10) 



Here 7Y,rl are terms that connect right and left movers, 
while 7i c hirai are chiral terms with all fcrmions moving in 
the same direction. 

Consider the couplings in Hrl- We have u>n = 
W22 = (convention), W12 — W21 (from Hcrmiticity), 
and A12 = A21 (from R <-> L symmetry). Thus there 

are 8 independent couplings: Wy^i ■> ^22 ■> an< ^ ^vi ■ 
Note that there are no four-fermion Umklapp terms in 
our two-band system. 

In the specific lattice model, we expand the interac- 
tions Eq. ([2]) in terms of the continuum fields and find 
"bare" values of the couplings: 



y n = v Q=a - 



Vq=o 



V 2k F1 

2 



A12 — Vq=q - 

An = 2V2fe F1 , 

A22 = 2V2fe F2 , 

A12 = 21^/2 j 



V, 



tt/2 



12 



Vk F1 -k f 



V, 



tt/2 



- 2K-/2 



Cfi — C22 — Vq=o - 

/-<p t t Vk F1 —k 
°12 - V Q=0 7T~ 



U 

~2 ' 



On — C22 — %U 
C 12 — 2Vk F1 -k F2 



(11) 
(12) 

(13) 

(14) 
(15) 
(16) 

(17) 
(18) 
(19) 

(20) 

(21) 
(22) 



Here V Q = £^=-00 V(x - x')e^ x ~ x > = V- Q , since 
V(x — x') = V(x' — x). We have also used explicitly 
k F i + k F 2 = -tt/2. 

The terms 7i c hirai renormalize "velocities" of various 
modes. In the weak coupling RG analysis, they only gen- 
erate higher order contributions and are therefore not im- 
portant. The RG equations below contain only couplings 
from H.RL- On the other hand, the chiral interactions 
are important in the intermediate coupling analysis to be 



done in Sec. IIII1 which is why we have listed their values 
as well. The on-site coupling U = V(x — x' = 0) appears 

explicitly in Cf/' 7 and because of our more careful 

treatment of the on-site interaction, which we first write 
as Uni (x)n±(x) and then insert the continuum fields (and 
bosonize in Sec. IIIII). 



B. Weak coupling Renormalization Group 



The RG equations in the two-band system are 



A 1X = - 
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(23) 
(24) 
(25) 
. (26) 
, (27) 
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(29) 
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Here O = dO / d£, where I is logarithm of the length scale. 
We have also defined 



A p/c 



X 



pi a 
11 



A 



p/cr 
22 



2A 



p/a 
12 



2lTVl 27TW2 n(vi + V2) 



(31) 



Details of our system enter through the band velocities 
Vi,V2, and the initial conditions Eqs. (fTTj) - (fl"8]) . 



C. Fixed point for stable C2S2 phase 

We are primarily interested in the stability of the two- 
band metallic phase with two gapless charge and two 
gapless spin modes - "C2S2" in the notation of Ref. l2ll . 
In the RG, this phase is characterized as having no diver- 
gent couplings. Before proceeding with detailed numeri- 
cal studies of the flow Eqs. (|23[) - (|3"01) . we can describe such 
stable C2S2 fixed point qualitatively: The charge sec- 
tor couplings reach some fixed values, A^, A22, A^, and 
are strictly marginal; they also need to satisfy A p * > 
(see below). The spin sector couplings approach zero 
from positive values, Af* = A22 = Af| = + , and arc 
marginally irrelevant. Finally, the "charge-spin" cou- 
plings W12 go to zero, = w i2 = 0, and are irrelevant, 
which is insured by the condition A p * > 0. Indeed, con- 
sider small deviations of comparable magnitudes for all 
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couplings and allowing only positive \° b . Since we have 

finite A p * > 0, first the vf^ will renormalize quickly to 
zero, without affecting significantly the other couplings. 
Then the \ a ah will renormalize to zero via slow marginal 
flows. 



D. Numerical studies of the flows 

We can solve the RG equations numerically for given 
initial conditions and check whether the couplings flow 
into the domain of attraction of the C2S2 fixed point or 
not. We use Mathematica to solve the flows up to long 
"time" I when the ultimate trends become apparent. 

If the couplings always remain of the same order as 
their initial values or approach zero, we say the cou- 
plings arc marginal or irrelevant and identify this as the 
C2S2 phase. The eventual trends here were discussed in 
Sec.liTCl 

On the other hand, if the magnitudes of some couplings 
grow significantly compared to the initial values, we say 
that the couplings are relevant and the C2S2 phase is 
destroyed. Thus, if either \\ x or AJ2 coupling becomes 
negative while w1 2 and w p 2 remain of the same sign, this 
A CT then runs away to large negative values and also in- 
duces the other couplings to diverge. Bosonizing the four- 
fermion interaction a 14 i 21 ' 25 (cf. Sec. we can see that 
two spin modes and one charge mode become gapped and 
we obtain so-called "CI SO" phase. The overall charge 
propagation mode remains gapless and the system is con- 
ducting. Note that we do not distinguish which coupling 
diverges faster in the formal flow Eqs. (j2"3"|) - (f3TJ)) . As dis- 
cussed in Ref. [2ll in the U — > 0+ limit one can separate a 
so-called "C2S1" case where one of the spin couplings di- 
verges qualitatively faster (but all couplings still diverge 
at the same £). We do not make such subtle distinction 
and call any runaway flow situation as C1S0 - all we 
want to know is that the two-band metal C2S2 became 
unstable. 

The RG flows are qualitatively similar for different 
points in the same phase, so we only show one repre- 
sentative picture for each case. Fig. [2] shows the flows 
in the C2S2 phase. The scale parameter £ is the x-axis, 
while logarithm of the couplings is the y-axis. In this 
way, we clearly see that the couplings separate into three 
groups, which is well explained by the C2S2 fixed point 
in Sec. Ill Cl the u>x2 now to exponentially rapidly, the 
A^ b flow to marginally slowly, while the \ p ab saturate. 

Fig. [3] illustrates the flows in the C1S0 phase. Here we 
use real values of the coupling as the y-axis and only show 
selected couplings, A^ 1; A£ 2 , w p 2 i an d w% 2 - We clearly 
see that these couplings diverge (and so do the other 
couplings not shown in the figure). 




FIG. 2: (Color online) Example of RG flows in the C2S2 
phase. The model potential is Eq. (|32[) with k = 1/2, 7 = 2/5; 
the band parameter is tz/ti = 0.9. We choose logarithm of the 
couplings to be the t/-axis and RG "time" i to be the x-axis. 
We see that w p ^ 2 " flow toward rapidly (irrelevant couplings); 
\ p ab saturate very fast (strictly marginal couplings) ; while \°t 
flow to slowly (marginally irrelevant). More generally, if we 
fix these k and 7 values, for ta/ii < 0.99 the flows are similar 
to those shown here and the phase is C2S2. 



Coupling 




FIG. 3: (Color online) Example of RG flows of selected cou- 
plings in the C1S0 phase. The model is the same as in Fig. [5] 
but with t%jt\ = 1.05. We see that the selected couplings 
diverge after some time. For example, once the A u and A22 
become negative while w p 2 and wf 2 remain positive, the RG 
equations (|23[) - (|30|) drive the An and A22 to -00 and in turn 
w p 2 and W12 to +00, and then all couplings diverge. More 
generally, if we fix 7 = 2/5, for ia/ii > 0.99 the flows are sim- 
ilar to those shown here and we call this C1S0 phase. Varying 
7, we obtain the phase diagram Fig. [4] 



E. Examples of phase diagrams with C2S2 metal 
stabilized by extended interactions 

For illustration in our paper, we consider the following 
interaction potential, 

v ^-^)={ KUe u ^i ; !*:£!>!} (32) 
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FIG. 4: Stabilization of the C2S2 metal by extended inter- 
actions. The model potential is Eq. (|32|l with k = 0.5. The 
non-interacting problem has one band for tiftx < 0.5 and 
two bands for t^/ti > 0.5, cf. Fig. [1] and we focus on the 
latter region. The limit 7 — > 00 corresponds to the Hubbard 
model with on-site repulsion only, and the C2S2 phase is sta- 
ble only over a narrow window t 2 /ti G [0.5 . . . 0.57].-i&2± The 
C2S2 region becomes progressively wider as we increase the 
interaction range I/7. 



Here U is the overall energy scale and also the on-site 
repulsion. The relative magnitude of the extended re- 
pulsion is set by some factor n < 1. Beyond one lattice 
spacing, the potential decreases exponentially with decay 
rate 7. For 7 — - > 00 we obtain the Hubbard model with 
on-site interaction only, while for small 7 the interaction 
extends over many lattice sites. 

We also consider the above potential but truncated at 
the 4-th neighbor. This tests robustness of our conclu- 
sions to modifications where the interactions have finite 
but still somewhat extended range, as may be preferable 
in numerical studies of such electronic models. 



1. Weak coupling phase diagram for potential Eq. 

The extended repulsion, Eq. (|32p . is in Fourier space 

Ksinh(7) 



Vq = U 



1-K 



cosh(7) — cos(Q) 



(33) 



For given model parameters, we use Eqs. (|TT]) - (|T5|) to 
set initial conditions. We follow the RG flows and iden- 
tify the phases as described above, thus mapping out 
the "weak coupling phase diagram" . Here and in the 
rest of the paper, we take k = 0.5. This is loosely 
motivated by the recent ab initio calculation^ for the 
k-(ET)2Cu2(CN)3 which gives the ratio of the nearest 
neighbor repulsion V\ = V(| x — x' \ = 1) to the on- 
site Hubbard U as V\/U ~ 0.43, while in our model 
V\/U = Ke~ 7 . The corresponding phase diagram show- 
ing stable C2S2 region is in Fig. 0J 



FIG. 5: Same as Fig. g] but for the potential Eq. trun- 
cated at the 4-th neighbor. 



We see that the C2S2 region becomes wider upon in- 
creasing the interaction range I/7. We can understand 
this qualitatively as follows. For fixed band parameters, 
when 7^0 the values of Vq for all non-zero Q approach 
[7(1 — k), while Vq—q ~ 2/cC//7 continues to increase. 
The corresponding contribution to A p is 



SA p = 



Vq=o 
2tt 

Vq=o 
2tt 



1 1 

V\ V 2 V\ 

(vi - v 2 ) 2 



V2 



ViV 2 (vi + v 2 ) 



(34) 



which is positive for any V\ 7^ v 2 and grows with increas- 
ing Vq= . Note also from Eqs. lfTT ]) -(fI% P that the V Q=0 
enters only in the X p ab couplings. Large bare value of A p 



makes the 



pi G 

11"' 



flows strongly irrelevant. Their effect on 



the \° ah flows is rapidly decreasing and expires. The \ a ah 
couplings start repulsive and stay so and eventually flow 
to zero via marginal flows. This argument is strictly true 
in the small 7 limit, while for finite 7 the interplay of 
different flows is more complex and requires numerical 
study as done in Fig. [¥] 



Weak coupling phase diagram for potential Eq. 
truncated at the J^-th neighbor 



Here, we truncate the interaction at the 4-th neighbor, 
so the Fourier transform is, 



Vq = U 



1 + 2k^V" 7 cos(nQ) 



The phase diagram in the weak coupling RG approach is 
shown in Fig. [5] We see that unlike the case without the 
truncation, the C1S0 phase opens again as 7 — > 0. Since 
we only include up to the 4-th neighbor interaction, Vq—q 
does not dominate over Vq^q even in the 7-^0 limit. 



6 



For k = 0.5 and 7 = 0, there is significant structure in 
Vq including sign changes as a function of Q, which can 
make bare spin couplings \°[ a ~ V2k Fa to be marginally 
relevant. Nevertheless, for intermediate 7 there is still a 
wide window of the C2S2 phase. 



III. WEAK TO INTERMEDIATE COUPLING: 
PHASES OUT OF C2S2 UPON INCREASING 
INTERACTION 

A. Harmonic description of the C2S2 phase 

Let us begin with a harmonic description of the C2S2 
metal. Technical steps and many details of the bosoniza- 
tion essentially follow Ref. and references therein. We 
write 



and "even" and "odd" flavor combinations, 



(37) 



with /! — p,cr. Similar definitions hold for the ip fields. 

We can now bosonize all four-fermion interactions 
Eqs. ([?])- (|T0j), First consider the spin sector. The C a ah 
terms give velocity renormalizations, while the \ G ah terms 
are written out in Sec. IVA of Ref. [IJ and are not re- 
peated here. We assume that the \" ah are marginally ir- 
relevant in the C2S2 phase. The fixed point Lagrangian 
has effectively decoupled boson fields B\ a and 82a with 
Luttinger parameters gi a = <?2ct = 1, dictated by SU(2) 
spin rotation invariance. 

The Lagrangian in the charge sector is 



(35) 



where tp and 8 are canonically conjugate boson fields and 
r\ are Klein factors. 

We define "charge" and "spin" boson fields, 
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lap/a — ~(^a| ± #aj.) , 



(36) 



C» = [d x & T - A-d x ® + d x * T -B-dx®] 

+-d x & T ■ a r * , 

7T 

where we defined @ T = (6 p+ ,8 p -) and <& T = ((p p+ ,(p p -). 
Matrix elements of A and B are: 
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,(39) 
(40) 
(41) 
(42) 
(43) 
(44) 



where 



Vl + V 2 



Vi - v 2 



(45) 



interaction^ 2 - 



J2V(x - x')n(x)n(x') V Q=0 f[p(x)f 



(46) 



The couplings \ p ab of the right-left mixing interactions 
H.ft L enter with opposite signs in A and B and directly 
affect Luttinger parameters, while the couplings G p ah of 
^chirai en ter with the same sign and give velocity renor- 
malizations. 

From the final expressions in terms of Vq , we see that 
the Q = component enters only in An. This can 
be understood by considering the Q = part of the 



where p(x) = 2d x 8 p+ /ir is the coarse-grained electron 
density. 

Note also that the —U/(2w) in the diagonal matrix 
elements is due to our more careful treatment of the on- 
site repulsion, which we first write as Un^{x)ni(x) and 
then bosonize. 

We obtain harmonic description of the C2S2 phase 
by combining the spin and charge sectors. The latter 
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two-mode system Cf has nontrivial Luttinger parame- 
ters, which can be determined from the matrices A and 
B (cf . Appendix |Aj . The fixed-point matrix elements 
will differ somewhat from the bare values above, but we 
ignore this in our crude analysis of the intermediate cou- 
pling regime. 

To complete the bosonization of the four-fcrmion in- 



teractions, Eqs. ([7|)-([T0]). the 



pi a 

'hi 



terms givi 



.14,25 



w 



{w p l2 J R i 2 J L i 



12 



JR12 ■ JL12 ) + H.C. 



(47) 



cos(2iy9p_)< 4u^ 2 cos(2i J 9 ._) — f cos(2# CT _) 



w a l2 cos(2<p CT _ ) + T cos(26» (T _ ) + 2f cos(20 CT ) } , (48) 



where Y = ^i| ? 7ij. 7 ?2T 7 ?2i- We see that W couples the 
charge and spin sectors. In the C2S2 theory described 
above, its scaling dimension is, 



A[W] = A[cos(2^ p _)] + I 



(49) 



where A[cos(2ip p -)] is evaluated in the Lagrangian C p , 
while the contribution 1 comes from the spin sector. For 
the C2S2 theory to be consistent, the W term must be 
irrelevant, A[W] > 2. Once the W renormalizes to zero, 
the charge and spin sectors decouple. We thus have pre- 
cise parallel with the weak coupling analysis of the C2S2 
fixed point in Sec. [TTJ 

On the other hand, if A[W] < 2, the W term becomes 
relevant and the C2S2 state is unstable. In this case, 
tpp— will get pinned and also the spin sector will become 
gapped. Only the mode remains gapless and the 

system is some CI SO conducting phase. 



B. Mott insulator driven by Umklapp interaction. 
Intermediate coupling procedure out of the C2S2 

The weak coupling analysis in Sec. |TT] misses the pos- 
sibility of gapping out the overall charge mode 6? p + since 
there are no four-fcrmion Umklapp terms allowed in the 
two-band system. However, the half-filled electronic sys- 
tem does become a Mott insulator for sufficiently strong 
repulsion. In the theoretical description, this is achieved 
by an eight-fermion Umklapp terir^ 



#8 = M4ii^4iii c R2^ c R2i c m c Lii c m c m + H - c -) 

= 2v 8 cos(46> p+ ) . (50) 

At weak coupling, this term has scaling dimension 
A[i?g] = 4 and is strongly irrelevant. However, from 
Eq. (|4l)|) we see that overall repulsive interaction stiffens 
the 9 P + mode and lowers the scaling dimension of H$. 
For sufficiently strong repulsion, A[_ff§] drops below 2 
and the Umklapp becomes relevant; 6 P + gets pinned and 
we obtain a Mott insulator. 



Our intermediate coupling procedure is as follows. Us- 
ing the harmonic theory of the C2S2 phase, we calculate 
the scaling dimensions A[W], Eq. (|4*9")h and A[i? 8 ] = 



A[cos(46* p+ )] from the Lagrangian C p , Eq. (|38[) . Details 
are described in Appendix [A] and calculations are done 
numerically in the end. 

If both A[W] and A[H S ] are larger than 2, the C2S2 
metal is stable. As interactions increase, eventually ei- 
ther W or H% becomes relevant. In general, there are 
two cases: 

1) If Hg becomes relevant first, we pin 9 p+ and en- 
ter "CIS2" Mott insulator. To be more precise, we can 
further qualify the label as "Cl[p— ]S2"; the remaining 
"charge" mode "p— " represents local current loop fluc- 
tuations and does not conduct. This is the spin liquid 
phase called Spin Bose-metal in Ref. and described in 
detail there. Exploring conditions for finding such phase 
in electronic models is our main goal here. 

2) On the other hand, if the W term becomes relevant 
first, we enter CISO conducting state with a spin gap 
(more precisely, "CI[p+]S0"). 

Some reservations arc in order. First, we use bare val- 
ues of the couplings in the A and B matrices, which is 
not accurate since the couplings experience initial flows, 
cf. Sec. Ill B[ Second, we consider only instabilities driven 
by changes in the harmonic C p theory as they translate to 
scaling dimensions of the H% and W terms, i.e., we effec- 
tively treat the latter as small. We also assume that the 
spin sector is near the fixed point with all marginally 
irrelevant and small. We will address these reservations 
after presenting results of the above procedure. Keeping 
these remarks in mind, we now describe how we analyze 
phases out of the Cl[p— ]S2 and Cl[p+]S0 in the same 
procedure. 



1. Instability out of Cl[p—]S2 driven by spin-charge 
coupling W 

In the present analysis focusing on the and "p— " 

fields, we can also crudely estimate the extent of the C1S2 
or C1S0 phases once either happens out of the C2S2. 

Suppose the Umklapp Hg, is relevant first and we are in 
the C1S2 phase. We still need to remember the W term 
since it can become relevant if we continue increasing the 
interaction strength. To estimate the scaling dimension 
of the W term, we assume now that the 8 p+ field is mas- 
sive and integrate out 9 p+ and ip p + . Mathematically this 
amounts to sending An — ► oo, and we obtain 



A[W; 



is pinned] = 



BnB 2 



Bf 2 



1. (51) 



This assumption is approximate but reasonable, since 
once the parameters are such that the system is in the 
CIS2 phase, the relevant Hg will grow and quickly stiffen 
the An m positive feedback loop. 

The CIS2 phase is stable if A[W] > 2, and this anal- 
ysis is similar to the stability analysis of the SBM in 
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Ref. 0- If A[W] drops below 2, the W term becomes 
relevant and the ip p _ field will be pinned, together with 
gapping out the spin sector, cf. Eq. (j48|) . The final result 
is some "COSO" phase, whose precise character depends 
on the details of the couplings w^" . This is studied in 
Sec. IVB of Ref. For the present repulsive electron 
model, we have w^ 2 ,Wi 2 > 0, so the resulting COSO is 
likely a pcriod-2 Valence Bond Solid (VBS)^ This con- 
nects to dimerized phase in the J\ — J2 spin chain ap- 
propriate in the strong interaction limit of the electron 
system. 



2. Instability out of Cl[p+]S0 driven by Umklapp H& 




0.5 0.6 0.8 1 1.2 1.4 1.6 



Suppose now the W interaction becomes relevant first. 
From Eq. (|4"8"|) . it is natural that tp p _ is pinned, the spin 
sector gets gapped, and we are in C1S0 phase. Here 
we postulate mass for ip p - (essentially sending B22 — > 
00) and calculate the effective scaling dimension of the 
Umklapp term, 



A [i? 8 ; <p p _ is pinned] = 4 



B11A22 



A n A 2 



A-12 



(52) 



If A[H 8 ] > 2, the C1S0 is stable. Once A[H$] drops 
below 2, the overall charge mode 8 P + is pinned and we 
obtain fully gapped Mott insulator COSO, which is likely 
the same period-2 VBS discussed earlier. 



C. Numerical results 

We consider the same models with extended density- 
density interactions as in the weak coupling analysis in 
Sec. Ill El parking ourselves initially in the C2S2 phase 
in Fig. U and Fig. [5] From the preceding discussion, we 
can obtain two phases out of the C2S2 upon increasing 
interaction strength - either Cl[p+]S0 or Cl[p— ]S2. To 
visualize the results, we imagine adding the overall in- 
teraction strength V as the z-axis to Fig. 0] and Fig. [5] 
We then project down which phase happens first for each 
such vertical line out of C2S2. Calculations are done nu- 
merically and the results arc shown in Fig. [5] and Fig. [9] 
In Fig. [7] we take a cut through Fig. [5] at 7 = 0.4 and 
show details of the phase diagram in the t2/ti — V plane. 



1. Intermediate coupling phase diagram for model with 
potential Eq. HIE) 

Fig. [5] shows results for the model potential Eq. (f3"2l . 
We can see that in two regimes 7 > 1.2 and 7 < 0.4 we 
exit from the C2S2 into the C1S2. The two limits can be 
understood analytically. 

In the large 7 case, we can replace all Vq by simply U. 



FIG. 6: Projection of phases obtained out of the C2S2 of 
Fig. [4] as we increase overall repulsion strength V, which we 
imagine to be the z-axis perpendicular to the page (Fig. [7] 
gives one cut at 7 = 0.4 with such V axis shown explicitly). 
The results are obtained in the intermediate coupling proce- 
dure as explained in the text. White region is C1S0 at weak 
coupling, cf. Fig. [4] and is not considered here. 



The matrices A and B defined in Eq. (|5g j) -([4"4" |) become 

(53) 



v v r 



We see that U only contributes to An. This monoton- 
ically "stiffens" the 9 P + (lowering A[ifg]) but "softens" 
the ip p - (increasing A [W]). Therefore we only expect the 
C1S2 phase out of the C2S2 as found in the numerical 
calculations. 

On the other hand, for small 7 we can see from Eq. (|33|) 
that Vq—o will dominate over Vq^q. Keeping only Vq = q, 
the matrices A and B become 



41', 



B 



V V r 

V r V 



(54) 



Thus the small 7 case has similar mathematical struc- 
ture to the large 7 case. The physical difference is that 
here the transition to the C1S2 is driven by the Vq = q 
instead of the on-site Hubbard U. Note also that since 
Vq—q ~ 2kU/j for 7 « 1, the transition requires only 
small values of U, which is why we can ignore all Vq^q 
compared to the band velocities. 

Now we consider a cut at 7 = 0.4 to see more details 
in the tilt\ — V plane. The results are shown in Fig. [7] 
Compared with the two limits 7> 1 and 7 C 1 above, 
all possibilities that we discussed out of the C2S2 are re- 
alized here. The C1S0 phase appears for tijt\ < 0.65 for 
some quantitative reasons. Various Vq are all of the same 
order, unlike the 7 <§C 1 case. At the same time, they have 
some non-trivial Q-dependence, unlike the 7 > 1 case, 
which is somehow enough to make the W term become 
relevant and preempt the Umklapp term. Note that for 
small interactions the scaling dimension of the W term 
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FIG. 7: (Color online) Intermediate coupling analysis of the 
model with potential Eq. ()32f) for k = 0.5 and 7 = 0.4. Here 
the horizontal range is equal to the extent of the C2S2 phase in 
the weak coupling analysis from Fig. [4] We start in the C2S2 
at small U . The boundary where the charge-spin coupling 
term W becomes relevant first is indicated with blue trian- 
gles and the system goes into the C1S0; the next stage where 
the C1S0 in turn becomes unstable and the system goes into 
the COSO is marked with green circles. The boundary where 
the Umklapp term H$ becomes relevant first is indicated with 
red squares and the system goes into the C1S2, which is the 
SBM phase of Ref. [l4l ; upon further increase of the interac- 
tion strength the C1S2 eventually becomes unstable and goes 
to the COSO at locations marked with black diamonds. Note 
that the discontinuity shown with dotted vertical line is not 
meaningful and is due to our crude analysis performed sepa- 
rately out of the C1S0 and C1S2; in either case, the final COSO 
is likely the same phase. Also note that the CI mode content 
is distinct in the Cl[p+]S0 (conducting) and CI[p— ]S2 (insu- 
lating) cases and any transition between them is first order. 
The C2S2 to C1S2 transition is Kosterlitz-Thouless-like. 



can be obtained from the weak coupling RG equations 
for the in Sec. Ill Bl by setting all \ a ah = (since we 
ignore the spin sector in the present procedure). Thus, 
A[W] = 2 + A p , where A p is defined in Eq. ([31]) ■ Since 
A p can only decrease under the weak coupling RG and 
the shaded C2S2 region in Fig. [4] was found to be stable, 
we expect A[W] here to increase with V for small V, in 
agreement with numerical calculations. However, we find 
that A[VF] eventually starts to decrease with increasing 
V and can become relevant before the Umklapp. This is 
a quantitative matter and comes from putting together 
all interactions Ti. p RL and W^ hiral , Eq. jT])-©, in the inter- 
mediate coupling procedure. Such numerical calculations 
give us that the C2S2 can exit into the C1S0 phase. For 
larger ^Ai > 0.65 in Fig. \7\ we obtain the sought for 
C1S2 spin liquid phase. 

This concludes the presentation of formal results 
within the particular procedure for intermediate scale 
analysis. Let us now think how to combine the weak 
and intermediate coupling approaches more realistically 
and see where our results are more robust. 



FIG. 8: Schematic merging of the weak and intermediate cou- 
pling results in the model regimes like in Fig. [7] in the whole 
range with £2/^1 > 0.5. In weak coupling, the C2S2 phase is 
unstable beyond the shaded region in Fig. [4] However, due to 
the crudeness of our intermediate coupling procedure, Fig. 
shows monotonic growth of the C2S2 phase with tz/ti past 
this instability. This discrepancy arises because our interme- 
diate coupling procedure completely ignores the spin sector. 
More realistically, we expect the C2S2 phase to peak some- 
where in the middle of the range shown in Fig. [7] and be 
bounded by the C1S0 for larger t%/t\. Similar considerations 
apply to the C1S2 phase, which is bounded by the COSO. 



First of all, in the weak coupling analysis the C2S2 
phase is unstable beyond the shaded regions in Figs. 2] 
and [5] However, this is lost in the specific intermedi- 
ate coupling procedure, which, when applied for small 
coupling, would give C2S2 essentially everywhere. For 
example, in Fig. [7] we see monotonic growth of the C2S2 
phase with t-xjt\ past the point where the weak coupling 
analysis predicts instability. The reason for this discrep- 
ancy is the complete neglect of the spin sector in the for- 
mal intermediate scale procedure. Indeed, in the weak 
coupling analysis, the instabilities manifest dramatically 
once one of the \ G aa becomes negative, causing runaway 
flows. This can happen even when the bare \ G aa are repul- 
sive because they are renormalized downwards and can be 
driven negative by the w^ 7 contributions in Eqs. (|26p ~ 
(|2T|) . where we assume w1 2 w1 2 > 0. Also, the A CT cou- 
plings feed back into the flow of Wy 2 a , so the RG flow 
behavior is even more complex. So far we have dealt 
with this inadequacy of the intermediate scale procedure 
by simply cutting it at the C2S2 boundaries determined 
from the weak coupling analysis. More realistically, we 
expect the extent of the C2S2 phase to peak somewhere 
in the middle of the range shown in Fig. [7] and decrease 
towards the right boundary. Similar considerations ap- 
ply to the C1S2 phase, which is likely confined within the 
same £2/^1 range as the C2S2. Therefore, the tijt\ — U lt\ 
phase diagram should be more like Fig. [5] 

We can also discuss our earlier reservation about using 
bare values of the couplings instead of some renormalized 
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values. Thinking about some RG treatment, we expect 
that crude patterns of how various couplings affect each 
other are likely similar at intermediate and weak cou- 
plings. Now if we formally take the flow equations from 
Sec. Ill Bl the outcome does not depend on the initial 
interaction scale, so we would conclude the C2S2 phase 
throughout the shaded region in Figs. 01 The weak cou- 
pling flow equations miss velocity renormalizations due 
to chiral interactions, but these are not expected to flow 
strongly and are treated reasonably in the intermediate 
coupling analysis. The fact that the couplings are now 
finite and comparable with bare band energies is also 
treated reasonably at intermediate coupling due to the 
power of bosonization, so the outlined forging of weak 
and intermediate scales seems appropriate. Finally, the 
Umklapp term that is missing in the weak coupling ap- 
proach will feed into stiffening of 6 p+ only, which is good 
for the first instability out of the C2S2 to be into the 
C1S2 spin liquid. 

We think that our conclusions are more robust for 
small 7 where the extent of the C2S2 phase is larger and 
also the longer-ranged potential is feeding precisely into 
stiffening the overall charge field 6 p+1 which is good for 
going to the C1S2 phase. On the other hand, results at 
medium to large 7 are likely less reliable, with different 
scenarios depending on quantitative issues. 



2. Intermediate coupling phase diagram for model with 
potential Eq. \32\l truncated at the 4~th neighbor 

Figure [5] shows results of the intermediate coupling 
analysis for the model with interactions truncated at the 
4-th neighbor, cf. Sec. Ill E 21 We have a rather simi- 
lar story to Fig. [SI except that the initial C2S2 region 
is bounded. Large part of the C2S2 phase exits into the 
C1S2 spin liquid upon increasing interactions, and our re- 
sults are probably more robust near 7 ~ 0.2 — 0.3 where 
the C2S2 has the largest extent along the t 2 /ti axis. 



IV. SUMMARY AND DISCUSSION 

To summarize, in this paper we consider electronic 
models for realizing Spin Bosc-metal (spin liquid) phase 
on the 2-leg triangular strip found in Ref. m spin- 1/2 
model with ring exchanges. We identify the SBM with 
the C1S2 Mott insulator of electrons. 

In Sec. HU we start with a two-band electron system, 
which is C2S2. Instead of considering only the on-site 
Hubbard-type repulsion ) 16 i 18 i 19 i 20 i 24 i 26 i 27 we study gen- 
erally longer ranged density-density repulsion. This is 
motivated in part by the expectation that real Coulomb 
interaction is not screened in Mott insulator materials, so 
further neighbor repulsion can be significant, as brought 
up by recent ab initio work 2 ^ for the spin liquid mate- 
rial k-(ET)2Cu2(CN)3. Using weak coupling RG analysis 
for the zigzag chain probleni ) 15 i 16 i 21 i 25 we find that such 
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FIG. 9: Same as Fig. [5] but for the model with interactions 
Eq. (|32|l truncated at the 4-th neighbor and starting out of 
the C2S2 of Fig. 



extended interactions open much wider window of the 
C2S2 metal compared with the Hubbard model. The 
main results arc shown in Fig. 3] and Fig. [5J In the first 
figure, we have essentially an independent control over 
the Q = part of the potential by allowing it to extend 
to far neighbors, and we identify the dominance of Vq = q 
as the main stabilizing force for the metal. In the second 
figure, we truncate interactions at the 4-th neighbor to 
check the robustness of our conclusions, in view that such 
models may be easier to explore using numerical DMRG. 
Our detailed quasi-lD considerations agree with the in- 
tuition that in real metals electronic pairing instabilities 
are suppressed by the long-ranged piece of the Coulomb 
interaction. Such widening of the C2S2 region by extend- 
ing the model interaction range is warranted if we want 
to bring the electronic ladder system closer to realistic 
situations in the 2D candidate spin liquid materials. 

In Sec. IIII| wc begin with stable C2S2 metal at weak 
coupling and use bosonization to extend the analysis to 
intermediate coupling by gradually increasing the over- 
all repulsion strength. Within effective bosonic theory, 
we identify potential instabilities of the C2S2 phase to 
spin-charge interaction W [Eq. (|4"5|) ] and Umklapp inter- 
action Hs [Eq. I|50p], The W can drive the system into 
Cl[/3-|-]S0 phase with spin gap but still conducting along 
the chain, while the Umklapp can produce Cl[p— ]S2 
Mott insulator with three gapless modes, which is the de- 
sired SBM phase. We calculate the scaling dimensions of 
the W and H$ terms in the harmonic theory of the C2S2 
metal using bare couplings in the charge sector and as- 
suming stability in the spin sector - this constitutes our 
naive intermediate coupling procedure. The calculation 
of scaling dimensions is described in Appendix and is 
done numerically in the end. 

We consider two cases depending on which of the terms 
W or H s becomes relevant first and apply similar in- 
termediate coupling approach inside the resulting phase. 
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Assuming strong field pinning by the already relevant 
term, we calculate the scaling dimension of the remaining 
term and estimate when it eventually drives the system 
into fully gapped COSO paramagnet (which is likely con- 
nected to the dimerized phase of the J\ — J2 Heisenberg 
model at strong coupling) . With the help of such admit- 
tedly crude analysis, we can map out the phase diagram 
in weak to intermediate coupling regime as illustrated 
in schematic Fig. [8] (based on more naive Fig. [7]). Fig- 
ures |5] and [5] summarize our results and show where the 
C2S2 metal goes to the C1S2 (SBM spin liquid) upon in- 
creasing overall repulsion strength. We conclude that the 
C1S2 phase is quite natural out of the wider C2S2 metal- 
lic region, in particular when driven by extended repul- 
sive interactions. It would be very interesting to confront 
our theoretical predictions with numerical DMRG stud- 
ies of such electronic models with extended repulsion. 

So far, we have approached the intermediate coupling 
Mott insulator from the weak coupling metallic side. One 
could try to attack the same problem starting from the 
strong coupling limit deep in the Mott insulator where 
Heisenberg spin-1/2 model is appropriate. As one nears 
the metallic phase, it becomes important to include mul- 
tiple spin exchanges in the effective spin Hamiltonian to 
better capture charge fluctuations in the underlying elec- 
tron systemic This is the motivation behind Ref. [3 
studying J\ — J2 chain with additional four-spin ring 
exchanges. The concept study Ref. [HI allowed arbi- 
trary variation of the ring coupling compared with the 
Heisenberg couplings. However, coming from an elec- 
tronic model these do not vary independently and more 
exchange terms are also generated. It would be inter- 
esting to pursue such approach systematically studying 
effective spin models with multi-spin exchanges for re- 
alistic electronic models to see if they harbor the SBM 
phase. We do not make such attempts here, but only give 
few simple observations on how the derivation of the spin 
model is modified in the presence of extended repulsion. 

First of all, for the two-spin exchanges, the familiar 
Hubbard model expression J rr > = 4t 2 , r ,/f7 is modified to 
J rr i = 4^ r ,/(Vo — V r - r '). The energy denominator is not 
simply the on-site U = Vq but also includes interaction 
potential between the two sites r and r'. For example, 
Ref. [H estimates Vl/V w 0.43 for the k-(ET) 2 Cu 2 (CN) 3 
spin liquid material, and this would significantly affect 
values of the exchange constants. Energy denominators 
for all virtual processes are similarly affected and take 
a form of a charging energy for the deviations from the 
background. Multi-spin exchange amplitudes are given 
by a product of electron tunneling amplitudes for a given 
virtual path divided by a product of such charging en- 
ergies in intermediate states along the path. Thus, the 
multi-spin exchanges may in fact be relatively more im- 
portant in systems with extended interactions. 

As an extreme example, imagine a very slow decrease 
of V(r — r') up to some distance R (and perhaps a faster 
drop thereafter). Then all exchange loops up to such 
radius R will have large amplitudes. The multi-spin 



exchanges encode the underlying kinetic energy of elec- 
trons, and our intuition is that this would like to retain 
some itinerancy in the spin degrees of freedom even when 
the charges are localized. From such strong to interme- 
diate coupling perspective, it appears that extended in- 
teractions would tend to stabilize the SBM spin liquid 
near the insulator-metal transition, similar to our con- 
clusion from the weak to intermediate coupling study in 
the quasi-lD models in this paper. It would be inter- 
esting to pursue such considerations more carefully and 
in realistic electronic models. We hope that our work 
will further stimulate numerical studies of such models 
on ladders and in two dimensions. 

Acknowledgments 

We would like to thank M. P. A. Fisher, I. Gonzalez, 
R. Melko, and D. N. Sheng for useful discussions and 
M. P. A. Fisher for stimulating this work and critical 
reading of the manuscript. This research is supported 
by the A. P. Sloan Foundation and the National Science 
Foundation through grant DMR-0907145. 

APPENDIX A: DERIVATION OF A[cos (40 p+ )] 
AND A[cos {2tp p -)\ IN C2S2 PHASE 

Equation (|38p gives quadratic Lagrangian for the 
charge sector. First, we redefine the fields which still 
satisfy the same commutation relations, 

= S-0i , * = S *i . (Al) 

Here S is an orthogonal 2x2 matrix diagonalizing the 
matrix A, 

S T -A-S=(^ l)^A D . (A2) 
The Lagrangian becomes, 

£p = h ^ dx&Ti ' Ad ' dx&i + dx ^ • sT • B • s • 8**1] 

+-d x &{ • d T 3> 1 . (A3) 
Define another set of conjugate fields, 

©i = ©2 , *i = v 7 ^ • *2 ■ (A4) 

V A -D 

We obtain, 

+ -d x &l ■ d T <f> 2 , (A5) 

where 

B' = Va^ ■ S T • B • S • . (A6) 
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We use the same trick to diagonalize matrix B': 



and scaling dimensions of the final fields. 



2 = R ©3 , * 2 = R # 3 , 

where R is an orthogonal matrix which satisfies, 

RT - B '- R =(o : B °i) sBi >- 

The Lagrangian becomes, 



(A7) 



(A8) 



C p = i- [d x ®% ■ d x & 3 + d x *% ■ B' D ■ d x $ 3 ] 

+-3 X &1 ■ 9 T * 3 • (A9) 

Now we can calculate the scaling dimension of cos (46> p +) 
and cos(2ip p -) from Eq. (|A9j) through relations, 



= S 



1 



R 3 , 



* = s v/a^ R * 3 , 



(A10) 
(All) 



A[e*** 3 ] 



fBL. 



A[e 1 * 3 ] 



4,/B 



(A12) 



D 



where the right hand sides mean corresponding diagonal 
matrix elements. Therefore, we find general form for the 
dimensions we are interested in, 



A[cos(40 p+ )] = 4^B[ 



aJb' 



j j S11R12 S12R 



22 



(A13) 



A[cos(2^ p _)] = 



/ AiS'2li?12 + VA2S22R22) 

7M 



{\U) 



where S„h and R n h are matrix elements of S and R. 
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